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ABSTRACT 

It has been claimed that Period-Luminosity relations derived from infrared observations of Large 
Magellanic Cloud (LMC) Cepheids are less dependent on the metallicity of the Cepheids. In this work, 
infrared observations of LMC Cepheids from the SAGE survey are combined with OGLE II optical 
observation to model and predict mass-loss rates. The mass-loss rates are fit to the data and are 
predicted to range from about 10~^^ to 10~''Mo/?/r; however, the rates depend on the assumed value 
of the dust-to-gas ratio. By comparing the relations derived from observations to the relations derived 
from predicted infrared stellar luminosities from the mass-loss model, it is shown mass loss affects 
the structure and scatter of the infrared Period-Luminosity relation. Mass loss produces shallower 
slopes of the infrared relations and a lower zero point. There is also evidence for non-linearity in 
the predicted Period-Luminosity relations, and it is argued that mass loss produces larger infrared 
excess at lower periods, which affects the slope and zero point, making the PL relations more linear 
in the wavelength range of 3.6 to 5.8 ^m. Because the dust-to-gas ratio is metallicity dependent and 
mass loss may have a metallicity dependence, infrared Period-Luminosity relations have additional 
uncertainty due to metallicity. 

Subject headings: Cepheids — circumstellar matter — Magellanic Clouds — stars: mass loss 
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1. INTRODUCTION 

Cepheids are powerful standard candles because they 
follow a Period-Luminosity (PL) relation. This relation 
has been determined using Cepheids in the Large Mag- 
allanic Clo ud (LMC) in optical and near infrared bands 
(jLanev k, Stobic 1994). LMC Cepheids also provide in- 
sight into stellar astrophysics as they have lower metal- 
licity relative to Galactic Cepheids, which has an effect 
on the pulsation of Cepheids. 

Recently, infrared PL relations have been derived 
using Spitzer observ atio ns from the SAGE pro gram 
(jMeixner et aH l2006l} _bv iNgeow fc Kanbud (|2008l ) and 
bv iFreedman et aD ()2008[ ). These infrared PL relations 
are important for extragalactic studies, and this will be 
even more so when the James Web Space Telescope be- 
gins operation. The infrared PL relations are powerful 
tools because metallici ty does not contribute significantly 
(jPreedman et al.|[2008[ ) and because the pulsation ampli- 
tude decreases in the infrared. However, IRAS observa- 
tions have found infrared excesses in Galactic Cepheids 
IPcasy 1988). Intcrferometric observations have also de- 
tected the existence of circumstell ar envelopes around 
a number of Gala c tic C epheids (iKervella et al.l l2006t 
iMerand et all l2006l 120071) in the K-band. The obser- 
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vations of infrared excess imply that there may be an 
additional uncertainty in infrared PL relations, and the 
excesses may play a role in LMC Cepheids. 

One proposed mechanism for generating circumstellar 
shells and infrared excesses about Galactic Cepheids is a 
stellar wind similar to those generated from other evolved 
pulsating stars. Mass loss is believed to generate shells 
about asymptotic giant branch stars [where the mass- 
loss rates are related to pulsation and dust condensation 
in the atmospheres, in both the Milky Way and the LMC 
(|Mattsson et al.l [2008D ] . One would not expect dust to 
form in the atmospheres of Cepheids because the temper- 
atures are greater than 1500 K, and dust-driving is not a 
plausible mechanism to generate a stellar wind. However, 
it has been argued that Galactic Cepheids use pulsation 
and s h ocks generated by pulsatio n to eject mass (iWillsod 
11989( 1. iNeilson fc Lested ()2008al ) developed an analytic 
model to study the affect of pulsation and shocks, in com- 
bination with radiative-line driving, on mass-loss rates 
in Galactic Cepheids, predicting rates of the order 10"^*^ 
to 1Q~'^ Mq/uv. It is argued that at some large distance 
from the surface of the Cepheid, the wind cools enough 
that a small fraction of the gas condenses into a dust shell 
that produces an infrared excess. The predicted infrared 
excesses from the theoretical model are consistent with 
both intcrferometric observations and IRAS observations 
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(jPeasvl [19881) . The analytic model was also applied to 
theoretical models of Cepheids with the metallicities of 
the Small and Large Magellanic Clouds, and the Galaxy, 
and mass loss was found to be significant in Mag ellanic 
Cloud Cepheids as well (jNeilson fc Lester '2008bl. This 
result, however, has not been tested by observations. 

The dust shells that surround Cepheids are opti- 
cally thin because they form at large distances from 
the Cepheids. For instance, a Cepheid with T^ff = 
6000 K will have a condensation radius of Tc/R* — 
0.5(rc/Teff)~^/^ = 16, assuming a condensation temper- 
ature of 1500 K. At the distance where dust forms, the 
dust shells are optically thin, and do not contribute to 
the extinction of starlight. 

Mass loss may also be a solution to the Cepheid 
mass discrepancy problem. The Cepheid mass discrep- 
ancy is the difference between Cepheid mass estimates 
based on stellar evolution isochrones and those based 
on stellar pulsation calculations. The mass discrep- 
ancy IS about 10-20% in the Milky Way (jCaputo et alj 
I200'5| ). abo ut 17-25% in the LMC and about 20% i n 
the SMC (iBrocato et all [200l iKeller fc WoodI [200l . 
iKeller fc WoodI (|2006D determined that the mass dis- 
crepancy increases as the metallicity decreases, and 
ICaputo et al.l (|2005| ) argue th at the d iscrepancy may also 
be a function of mass, though iKelleri (,2008,) provides ev- 
idence against this result. For mass loss to be a solution 
to the mass discrepancy, a Cepheid that starts on its first 
crossing with mass of ^Mq would need to lose about IMq 
or have an average mass-loss rate in the range of 10~* 
to lO-^Mo/yr. 

The purpose of this work is to model infrared excess 
in LMC Cepheids using SAGE observations in the IRAC 
b ands combined with O GLE II observations of B,V, and 

I (jUdalski et al.lfl99 9a,bl. The next section outlines the 
observations and the model describing circumstellar dust 
created in a stellar wind that causes infrared excess. The 
process for determining the mass-loss rates is also de- 
scribed. The results are given in Section 3 and the pre- 
dicted infrared PL relations are described in Section 4. 
The fifth section will explore possible driving mechanisms 
for mass loss in LMC Cepheids, testin g if the mass loss 
behavio r is similar to that proposed in lNeilson fc Lester! 
(I2OOM . 

2. THE DATA AND MASS LOSS MODEL 

We use OGLE II and SAGE observations of LMC 
Cepheids to determine mass-loss rates. The OGLE II 
data are for B, V and I magnitudes while the SAGE mag- 
nitudes arc in IRAC bands at wavelengths 3.6, 4.5, 5.8 
and 8.0 /im. The SAGE data we used are adopted from 
iNgeow fc Kanbui] (|2008f l. which consist of 730 OGLE 

II LMC Cepheids with logP > 0.4. However we only 
use 488 of these Cepheids that have at least 3 IRAC 
bands and 2 of the BVI bands from this dataset. The 
SAGE data are compiled by matching the position of 
OGLE II Cepheids with positions of infrared sou rces in 
the SAGE observations. INgeow fc Kanburj ()2008[) match 
sources if they are within 3.5 arcseconds of the position 
of the OGLE II Cepheids. Figure [T] show the infrared 
PL relations constructed using the IRAC magnitudes. A 
number of (mostly short period) Cepheids appear to de- 
viate from the infrared PL relations, implying there is 
some infrared excess. 
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Fig. 1. — The observed brightnesses of LMC Cepheids in the 
infrared as a function of period. 

The are three possible causes of the infrared flux ex- 
cess: (1) blending of stars in SAGE observations, (2) 
false matches of the infrared sources to the OGLE II 
Cepheids, (3) or circumstellar dust shells forming at a 
significant distance from the Cepheids in a stellar wind. 
It is possible that false matches contaminate the sample, 
and we check this in Figure [2] where the magnitude resid- 
uals of the IR Peri o d-Lu minosity relations determined by 
INgeow fc Ka^Ib^ ^M) are shown as a function of the 
separation between the OGLE II an d SAGE positions. 
Although the search radius used in INgeow fc Kanburl 
([2OO8) is rather large, most of the matched objects have 
a separation less than 0.77 ar cseconds [as demonstrate d 
in the Figure 1 and Table 1 in INgeow fc KanbuH (|2008D ]. 
The potential for false matches is onl y significant for a 
small fraction of the total sample of iNgeow fc Kanburl 
(2008"), and in the sample used here, only ten of the 488 
Cepheid matches have a separation greater than 1.3 arc- 
seconds, which is the separation where the matches al- 
most all have large residuals to the fit of the PL relation. 
There is another asymmetry apparent in Figure ID where 
the residuals have a separation less than 0.4 arcseconds. 
In this work, we keep the two samples, but identify them 
in the figures when important. 

We can test if the infrared excess is due to mass loss 
by calculating the sum of the infrared luminosity of the 
Cepheid and the luminosity of the dust that is generated 
in the wind, given by 

3 <a'>Qt Md 
47r < a-^ > p Vd 
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Fig. 2. — The apparent magnitude residuals of the fits of the PL relations from lNgeow h Kanburl 12008) as a function of the separation 
of the positions of OGLE II Cepheids and the infrared sources for (Left) 3.6 and 4.5 /im and (Right) 5.8 and 8.0 /im. 



X / B,[Td)[l-W{r)]dr. (1) 

JR., 

The dust shell luminosity is proportional to the ratio of 
the mean cross section, < >, and volume, < >, 
of the dust particles, and inversely proportional to the 
mass density of dust particles, p. The dust mass-loss 
rate is given by Md^ the dust velocity is Vd, and 
is the absorption efficiency. The term in the integral 
is the product of the blackbody radiation of the dust 
with temperature and the geometric dilution factor 
W{r) = [1 - - (i?,/r)2]/2 at a distance r from the 
surface of the Cepheid. 

To predict the dust shell luminosity, we need to spec- 
ify the properties of the dust. The dust grain size is as- 
sumed to range from 0.005 /xm to 0.25 ^m, which yields 
a value < a? > / < >w 40 pm based on the method 
of iMathis et al.l ()1977[ ). Assuming the dust is primarily 
graphite, the absorption will be concentrated at optical 
wavelengths, giving an absorption efficiency of « 2. 
This further implies that the mean density of the grains 
is p — 2.2 g/cm?. The dust velocity, equivalent to the 
terminal velocity of a wind, is about 100 km/s, approx- 
imately equivalent to the escape velocity of a Cepheid. 
The integral is computed from the surface of a Cepheid 
but dust does not form in the wind until the material is 
at a condensation distance Tc — (i?*/2)(T*/1500iir)^/^, 
where dust condenses at a temperature of 1500K. The 
dust temperature at distance r from the star, greater 
than the condensation distance, is Td{r) — T,VK(r)^/^. 
This leaves the dust mass-loss rate, stellar radius and 
effective temperature as unknowns in Equation [T] 

The gas mass-loss rate is found by assuming a dust-to- 
gas ratio. The typical ratio assumed for the Milky Way 
ISM is 1/100, and this ratio was use d in previous stud- 
ies for mass-loss in Galactic Cepheids ([McAlarv fc WelchI 
[1986) . The dust-to-gas ratio in the LMC must be signif- 
icantly smaller than 1/100 because the formation of dust 
depends on the metallicity of the gas. For this work, 
a value of the LMC dust-to-gas ratio is assumed to be 



1/250 found by scaling the Milky Way dust-to-gas ratio 
by the ratio of the average LMC metallicity oi Z = 0.008 
to the standard solar metallicity Zq — 0.02. This choice 
of dust-to-gas ratio leads to a gas mass-loss rate that is a 
lower limit. The dust-to-gas ratio in the LMC has been 
observ ed to be approxim ately one quarter the Galactic 
value (jClavton fc Martin ] |1985D to about on e tenth the 
Galactic value (jWeingartner fc Drainell2001[ ). Therefore 
a gas mass-loss rate may be smaller than would be pre- 
dicted using other dust-to-gas ratios. The dust velocity, 
which is also the terminal velocity of the gas wind, is cho- 
sen to be approximately the escape velocity, but the dust 
velocity may range from about 75-150 km/ s, leading to 
an uncertainty of about 50%. 

The mean luminosity of a Cepheid at frequency v is 
given by 

,Star 

meaning the stellar luminosity is dependent on the ra- 
dius and effective temperature. We check if a blackbody 
is a reasonable approximation by comparing blackbody 
B, V, I brightnesses wit h the B,y, and I from a Atlas 
stellar atmosphere mode lKurucd (|1979t ) with an effective 
temperature of 6000if and \ogg = 1. The model at- 
mosphere is consistent with Cepheid properties and the 
B,V, I agree with blackbody estimates to within a few 
tenths of a magnitude. The infrared observations are well 
approximated by a blackbody brightness because these 
wavelengths are in the tail of the blackbody function at 
6000 K. The different predicted brightnesses will affect 
the uncertainty of the mass-loss model but not greatly. 
This leaves three unknown variables for fitting the ob- 
servations: the dust mass-loss rate, the stellar radius 
and the effective temperature. The effective tempera- 
tu re may be de termined using the relation determined 
bv lBeaulieu"er al. (2001) 

logToff = 3.930122 -I- 0.006776 logP - 0.2487(y - /)o. 

(3) 

The relation is dependent on the pulsation period and 
de-reddened color. 
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To compare the predicted total luminosity (stellar plus 
shell) of the LMC Cepheids with the observed fluxes, we 
need to adopt a distance modulus. Our choice is 18.5±0.1 
for the L MC, which is approximately the mean distance 
modulus (ICatelan fc CortesI 120081 : IClement et aD l2008l . 
for example), although the value ranges from 18.4 to 18.7. 

The errors for the optical and infrared observations are 
likely to be negligible compared to other uncertainties, so 
we assume that the error for each optical wavelength is 
0.1 magnitudes, due to the uncertainty of the distance 
to the Cepheids to determine the absolute magnitudes 
and the error for the infrared observations is 0.2 mag- 
nitudes based on the distance uncertainty. The main 
sourc es of error are the thickness of the LMC (|Lah et al.l 
[2005'), error in the observations themselves and that the 
infrared observations may not be the mean brightness of 
Cepheids. In fact, a thickness of about 5% of the distance 
to the LMC corresponds to an error of 0.1 mag. 

The BVI observations are also de-reddened while it 
is assumed that the extinction of the infrared light is 
negligible. The color excess E{B — V) is given in the 
OGLE II data for each Ccpheid and is on average « 0.15. 
The extinction for the BVI is calculated in the same way 
as in lUdalski et al.l (|1999bD , while the extinction in the 
infrared is significantly less than 0.1 magnitudes. 

3. QUALITY OF FIT OF MASS-LOSS RATES 

Having described the method for predicting stellar 
fluxes over the seven bands, we fit the mean radius and 
dust mass-loss rate, and hence the gas mass-loss rate, 
using fitting, for each Cepheid in the sample. The 
value of is given as 

2 _ 1 (Theory) - M^, (Obs) \ ^ 

where N is the number of observations and / is the 
number of unknowns. The fits are calculated for a 
rang e of stellar radii ba sed on the Period-Radius rela- 
tion (" Gieren et aLlUQOgf ) 

log i?/i?0 = 0.68 logP(d) -f 1.146. (5) 

and varying that radius by ±20%. The dust mass-loss 
rate is fit for a specific mean radius by finding a minimum 
value of using a root-finding algorithm. Before we 
attempt to fit the dust mass-loss rate, we first apply the 
X^ fitting to the range of mean radii alone as a reference. 
The values of x^ and predicted radii in this method are 
shown in Figure [31 

Fitting the radius alone to the observations seems to 
provide a reasonable fit, but one may argue that the un- 
certainty of the fit is related to the fact that the SAGE 
observations are single epoch and are not mean bright- 
ness at these wavelengths. We test this by fitting 
the radius and a quantity dm. This variable represents 
the difference between the observed infrared magnitude 
and the mean brightness of the Cepheids. The values of 
dm are assumed to vary from —0.5 to 0.5 magnitudes 
dimmer, where a negative value of dm implies the mean 
brightness is brighter than what is observed. The cho- 
sen range of values are exaggerated as the full brightness 
amplitude of a Cepheid at IRAC wavelengths is < 0.4 
the amplitude in the visible. The largest full amplitude 
in the visible is 1.2 magnitudes, meaning the infrared 



variable dm will be —0.25 < dm < 0.25 in reality, and 
for the majority of Cepheids the range of dm is much 
smaller. We show in Figure [3| the predicted values of 
dm as a function of period. If the hypothesized infrared 
excess were due primarily to the fact that the infrared 
observations are single epoch then about 50% of the sam- 
ple would predict dm < and 50% > 0. The results in 
Figure [4| show a preference for dm > 0, with 68% of the 
total sample having dm > and 13% preferring dm = 0. 
Therefore, we take this as proof of an infrared excess in 
the SAGE sample. 

The next step is to determine the best-fit radius and 
mass-loss rate for the sample of Cepheids. The x^ fits 
are presented as a function of pulsation period in Fig- 
ure [5| along with the gas mass-loss rate. The mass-loss 
rates range from 10^^^ to 10~^MQ/yr, with the values 
of ranging from about 1.1 to about 23. The large ma- 
jority of Cepheids appear well fit by a mass-loss model 
that forms dust at some distance from the surface of the 
Cepheid. These predicted mass-loss rates are significant, 
10~^ MQ/yr, and, depending on the dust-to-gas ratio 
and dust grain properties, the gas mass-loss rates may 
be an order of magnitude larger or even more. These gas 
mass-loss rates are the minimum mass-loss rates for the 
LMC Cepheids. 

We also quantify the uncertainty of the mass-loss rates 
caused by the unknown phase of the infrared observa- 
tions by computing the best fits for the mass-loss 
rate, radius and the quantity dm that was defined ear- 
lier, and we compute the values of 5\n{M)/d{dm) for 
a random, with respect to period, subsample of 100 of 
the Cepheids. This error is a function of both the pul- 
sation amplitude, and the mass-loss rate. We show the 
values of 6\n{M)/6{dm) as a function of M in Figure [6l 
The uncertainty is related to the mass-loss rate, and, as 
one would expect, the uncertainty of the pulsation am- 
plitude is less important for larger predicted mass-loss 
rates. We highlight the boundary where the uncertainty 
of the mass-loss rate is 100% for a pulsation amplitude of 
0.5 magnitudes which represents the maximum infrared 
pulsation amplitude. This implies that Cepheids with 
predicted mass-loss rates > 10~^ MQ/yr have infrared 
excesses that cannot be explained solely by pulsation 
phase, and a number of Cepheids with mass-loss rates 
< 10~^ M^/yr have uncertainties S\n{M)/6{dm) < 4. 
As the IR pulsation amplitudes become known we will 
be able to probe smaller mass-loss rates with more cer- 
tainty. 

The model is also tested by comparing the two param- 
eter fits, with the mass-loss rate and radius as the two 
degrees of freedom, to fits using just the radius as the 
only degree of freedom. We use the F-test to quantify 
the significance of the mass-loss model . The F-test is de- 
scribe d by Kanbur & Ngcow (2004a); iNgeow fc Kanbuil 
()2008l and references therein) and for each Cepheid we 
calculate the value of F. For the majority of Cepheids, 
the fit of the radius has only six degrees of freedom while 
the mass-loss model has five. Some of the Cepheids have 
one less degree of freedom for each model respectively. 
The values of F are shown in Figure [7| against the val- 
ues of x^ for the mass-loss model. The values of F > 5 
mean that we can state with 95% confidence that the 
mass-loss model is significant relative to fitting the ra- 
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Fig. 3. — (Left) The fit of the observations w ith best— fit mean radius only. (Right) The predicted radii of the sample of Cepheids 
with the Period— Radius relation l IGieren et al.|[T999l) shown as a dotted line. The circled points are the Cepheids that are considered most 
likely to be false associations while the triangles represent the Cepheids with separation between OGLE II and SAGE coordinates of less 
than 0.4 arcseconds. 
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Fig. 4. — The best— fit values of dm, the predicted differences 
between the mean brightness and observed brightness at infrared 
wavelengths. Circled points are the likely false associations and 
triangles represent separations less 0.4 arcseconds. 

dius alone. This provides no information into the pos- 
sibihty of blending, for instance, so we also take a cut 
of models with a value of < 5. This corresponds to 
the upper left quadrant of Figure [7] which contains 44 
Cepheids. Therefore we state with 95% confidence that 
approximately 9% of the sample of LMC Cepheids have 
circumstellar dust shells caused by stellar winds. 

It is shown that 44 of the Cepheids are consistent with 
the mass-loss model, implying that the remainder of the 
sample is consistent with no mass loss. However the value 
of _F is a function of the mass-loss rate, with F increasing 
with M, as shown in Figure [71 Cepheids with values of 
F < 9 and < 5 in the lower left quadrant of Figure [7] 
with x^ < 5 are consistent with mass loss but the obser- 
vational errors are too large to state with certainty that 
all LMC Cepheids are undergoing mass loss. This sug- 
gests that the model needs to be tested with time series 
infrared observations to constrain the pulsation ampli- 
tude and reduce the uncertainty of the infrared observa- 
tions. In the next section, we test what effect mass loss 
might have on the infrared PL relation if the hypothesis 
is correct. 



4. THE EFFECT OF MASS LOSS ON INFRARED PERIOD 
LUMINOSITY RELATIONS 

It has been postulated that mass loss generates dust 
in a circumstellar shell surrounding a Cepheid and this, 
in turn, affects the infrared PL relation. By using the 
predicted stellar luminosities of the sample of Cepheids, 
we compute the stellar PL relations, m\ = a log P + b, 
and c omp are them with re sults from iNgeow fc Kanbuil 
(|2008D and'Fr eedman et all |2008). We also test the data 
for non-linearity in the infrared PL relations. In the fit, 
we do not include the Cepheids with large separations 
that were noted in the previous section, however, we do 
include those with separation less than 0.4 arcseconds be- 
cause they have randomly distributed mass-loss rates, in 
Figure [HI and hence have randomly distributed infrared 
excesses. 

The predicted stellar luminosities of the LMC Cepheids 
in the mass-loss model are shown in Figure |51 together 
with a comparison of the line ar and non-linear PL re - 
lati ons with the rela t ions f rom lNgeow fc Kanbuil ()2008D 
and lFreedman et al.l ()2008[ ). The non-linear PL relation 
is defined as two linear relations, the first for the period 
range of 1 t o 10 days while the second is for the longer pe- 
riod range (|Ngeow et all 120051) . The slopes, zero-points 
and dispersions of the fits for the linear and non-linear 
fits are given in Table [1] 

The linear PL relati ons have smaller d ispersion than 
the PL relations from iNgeow fc Kanburl (2008) for the 
two longer wavelengths and slightly larger dispersion for 
the two shorter infrared wavelengths. There are also no- 
table differences in the slopes and the zero points. The 
zero points of the predicted 3.6 and 4.5 PL rela- 
tio ns are approximately th e same as those determined 
bv INgeow fc Kanbu r (2008), while at longer wavelengths 
the differences are significant and also the predicted zero 
points ten d to be a little brighte r than the zero-points 
found by iFreedman eFall (|2008f) . This is consistent 
with mass loss causing larger infrared excesses at longer 
wavelengths. The predicted slopes range from —3.14 at 
shorter wavelengths to —3.18 at longer wavelengths. This 
is a small change of slope as a function of wavelength, and 
it is roughly consistent with a constant slope within the 
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Fig. 5. — (Left) The of the observation with best-fit mean stellar radius and mass-loss rate. (Right) The gas mass-loss rates of the 
sample of Cepheids, binned into different groups. Those points circled are the possible false associations. 
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Fig. 6. — The uncertainty of the mass-loss rates due to the dif- 
ference between the observed infrared brightness and the unknown 
mean brightness of the Cepheids as a function of mass-loss rate. 
The horizontal line refers to an uncertainty of 100% for a full pul- 
sation amplitude of 0.5 mag. 



TABLE 1 

Best Fit Parameters for Predicted PL Relations 



Type 


X {/J-rn) 


Slope 


Zero Point 


Dispersion 


Linear 


3.6 


-3.145 ±0.024 


15.993 ±0.017 


0.110 




4.5 


-3.159 ±0.023 


15.921 ±0.017 


0.108 




5.8 


-3.170 ±0.023 


15.924 ±0.017 


0.107 




8.0 


-3.181 ±0.023 


15.929 ±0.017 


0.105 


Non- 


3.6 


-3.248 ± 0.038 


16.057 ±0.025 


0.107 


Linear 


4.5 


-3.259 ±0.038 


15.983 ±0.025 


0.105 


P < lOd 


5.8 


-3.268 ±0.037 


15.984 ±0.025 


0.104 




8.0 


-3.276 ±0.037 


15.988 ±0.024 


0.103 


Non- 


3.6 


-2.971 ±0.123 


15.815 ±0.146 


0.125 


Linear 


4.5 


-2.989 ±0.121 


15.747 ±0.143 


0.122 


P > lOd 


5.8 


-3.005 ±0.119 


15.754 ±0.141 


0.120 




8.0 


-3.019 ±0.118 


15.763 ±0.140 


0.118 



errors given n Table[TJ The slopes from lNgeow fc Kanbuii 
([2008') show the opposite behavior with the slopes be- 
coming less st eep with longer wavelen gth, contrary to the 
arguments in iFreedman et al.l (|2008f) who maintain the 



slope of the PL relation should be steeper as a function 
of wavelength and approaches an asymptotic limit. In 
both cases the uncertainty of the slope is similar, rang- 
ing from 0.017 to 0.048 with increasing wavelength for 
iNgeow fc Kanbur' f2008), and about an average of 0.03 
for Frccdman et al. (2008). Therefore the slopes from 
these two works and those predicted here do not agree 
within the uncertainty. 

This shift in the behavior of the slopes of the infrared 
PL relations is due the removal of infrared excess caused 
by mass loss. This is best seen at 8.0 jum where the 
slope changes from about —3 (|Ngeow fc Kanbui1l2008f ) to 
—3.18 when the contribution of brightness due to mass 
loss is removed, although the observed slope of —3 may 
also be due to incompleteness of the data at the faint end 
of the PL relation. The implication is that mass loss in 
short period Cepheids contributes significant luminosity, 
increasing the zero point and causing a shallower slope 
because the majority of the Cepheids in the sample have 
periods less than ten days. The shallower slope implies 
that mass loss contributes fractionally less to the total 
infrared luminosity at longer periods because the stellar 
luminosity is already so large. 

The analysis in the previous section showed that the 
mass loss hypothesis is statistically unique from fitting 
only the radius of the Cepheid for about 44 Cepheids in 
the sample, or conversely that the majority of Cepheids 
in the sample are statistically consistent with zero mass 
loss. From this realization, we wish to compare the ob- 
served infrared brightness that we fit our model with and 
the predicted infrared brightness of the Cepheids. The 
comparison is shown in Figure [5] for the 8.0 /im data. At 
this wavelength the differences is most apparent because 
a dust shell contributes a larger fraction of the total flux 
at longer wavelengths. We also compute best-fit linear 
relations for the observed data where the Cepheids with 
large separation and the 44 Cepheids where mass loss is 
shown to be likely are not used in the fitting. This rela- 
tion is mg.o^m (Observed) = —2.905 log P + 15.530 with 
a standard deviation of 0.219. For comparison, we de- 
rive the best-fit data us ing all of the 8.0 /xm data from 
INgeow fc Ka^A^ ([2008h and find wg.o^m (Complete) = 
-2.473 logP-h 15.058 with a standard dev i ation of 0.602. 
The relations given in iNeeow fc Kanbuii (l2008h are de- 
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predicted mass-loss rates. The circled points represent those Cepheids that are likely false associations of the infrared observations and 
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termined using an iterative fitting method where a best- 
fit relation is determined and then any Cepheids with 
a brightness that is more than 3cr different are removed 
and a new PL relation is computed and the process re- 
peats until the PL relation converges. Here, we compute 
the PL relation using all of the data without the iter- 
ative approach, The linear relations are shown in Fig- 
ure [HI Although we are only able to confidently state 
that 44 of the Cepheids are consistent with the mass- 
loss model, we note that there are significant differences 
between the predicted stellar brightnesses and observed 
data, and that these differences are reflected in the in- 
frared PL relations. It is also interesting that these 44 
Cepheids are observed to be brighter than the majority 
of the sample but there are a number of Cepheids with 
similar brightness that are statistically consistent with 
zero mass loss. 

The data are tested for non-linearity in the infrared PL 
relations. The hypothesis that the infrared PL relations 
are non-linear is tested w i th the F-test, as describ ed in 
iKanbur fc Ngeowl (j2004a[ ): iNgeow fc Kanbuil (|2008l . and 
references therein) by comparing a PL relation of the 
form 

_(a\ogP + b logP<l 
\c\ogP + d logF>l 

with the standard linear PL relation with two degrees 
of freedom. If the value of F > 3, the PL relations are 
non-linear with 95% confidence. Our values of F, with 
increasing wavelength, are 7.89, 5.93, 5.78, and 5.50. The 
predicted stellar PL relations are thus consistent with 
being non-linear with a period break at 10 days. The 
predicted slopes and zero-points are given in Table [TJ 
However there are two possible sources of error. The 
first is that we are assuming blackbody radiation that 
ignores any infrared absorption lines that may affect the 
structure of the PL relations. The second is that there 
are significantly less data for periods greater than 10 days 
(approximately 50 data points). 

The non-linearity is related to the fact that the lumi- 
nosities are given by the effective temperature; using the 
OGLE II data to derive effective temperatures may cause 



a non-linear Period-Temperature relation because of 
non- hnearity in the OGLE II (T^ — /) Period-Color rela- 
tion (|Kanbur fc Ngeowll2004bf ). This non-linear Period- 
Temperature relation causes non-linearity in the infrared 
predictions. This implies that the PL relations given by 
only the stellar component is non-linear in the wave- 
length range of 3.6 to 8.0 /xm based on blackbody ar- 

fments, contradic ting the r e sults of iNgeow fc Kanburl 
)08) and Nge ow fc KanbuH (|2006D for the K-band PL 
relation. There are two plausible reason s why this con- 
tradiction is found. iKanbur et al.l ()2004D argue the non- 
linearity is due to the hydrogen ionization front (HIF) in- 
teracting with the photosphere, causing significant tem- 
perature variations in the layers of the Cepheids that 
emit mostly in the optical; at longer wavelengths this in- 
teraction becomes less significant. This implies that the 
mean effective temperature at shorter periods is affected 
by the HIF while at longer periods the effective temper- 
ature is just what would be expected for a non-pulsating 
star. At infrared wavelengths, most of the radiation is 
emitted higher in the stellar atmosphere farther from the 
effects of the HIF, which leads to a more linear PL rela- 
tion and is hence more dependent on the Period-Radius 
relation, which is linear. This would explain why the val- 
ues of F for the non-linear relations decrease with longer 
wavelength, the IR PL relations are becoming more con- 
sistent with a surface brightness related to the linear 
Period-Radius relation. The non-linearity of the pre- 
dicted data may just be reflecting the non-linearity in 
the optical wavelengths because the variations of the ef- 
fective temperature over the pulsation period is ignored. 

The second possibility is that mass loss causes larger 
infrared excess for shorter period Cepheids than for 
longer period Cepheids even though the mass-loss rates 
are similar for short (P < lOd) and long {P > lOd) pe- 
riod Cepheids. The short period Cepheids have smaller 
radii leading to smaller, more dense circumstellar shells. 
The more dense shells cause greater infrared excess. This 
greater infrared excess in shorter period Cepheids makes 
them appear brighter on average, which increases the 
zero point of the infrared PL relation; because the rel- 
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ative infrared excess decreases with longer period, the 
slope of the PL relation will appear shallower, with 
the effect being more prominent at longer wavelength. 
This idea may explain the marginal linearity found in 
the K-band PL relation. Infrared excess in Galactic 
Cep heids has been observed using K-band i nterferome- 
try (jKervella et al.ll2006t iMerand et aT]|2006l l2007i ) so it 
is likely infrared excess plays a role in the LMC Cepheids 
at this wavelength. This would imply that the K-band 
PL relation is actually non-linear and this non-linearity 
is being masked by the infrared excess. This argument 
also explains the results of the tests of non-l i nearit y in 
the IRAC PL relations in iNgeow fc Kanbuii (|2008f ). in 
particular the non-linear PL relation at 8.0 fim. The 
authors found that the slope of the non-linear PL rela- 
tion for P < lOd is shallower than the slope of the linear 
relation at 8.0 /im with a more luminous zero point. The 
non-linear PL relations for P < lOd in the optical and 
near-IR all display the opposite behavior with respect 
to the linear PL relations. This non-linear relation at 
8.0 /im is due to the same process that causes the other 
IRAC PL relations to appear linear except the process is 
more significant at longer wavelengths. 

It has been shown that mass loss provides a significant 
contribution to the infrared brightness of LMC Cepheids 
and affects the structure of infrared PL relations. With- 
out the contribution of mass loss, the slopes of the linear 
PL relations are steeper with increasing wavelength al- 
beit at a small rate diffe ring from the slop e s bec oming 
more shallow as found by I Ngeow fc Kanbun (|2008f ) . Ap- 
plying the F-test to the predicted data implies that the 



infrared PL relations are non-linear, though this result 
requires further testing. However most of the Cepheids 
have predicted mass-loss rates that are statistically con- 
sistent with zero implying this result is preliminary and 
needs to be tested further with more data with smaller 
uncertainties. 

5. WHAT IS THE DRIVING MECHANISM? 

Up to this point, we have investigated the ability of 
mass loss to match the OGLE II and SAGE observa- 
tions of LMC Cepheids, and how the resulting estimates 
of infrared excess affects the structure of the PL rela- 
tions. This has been done without assuming a driving 
mechanism of the Cepheid wind. There are a number of 
possible methods for stars to drive mass loss, but only 
two are likely for Cepheids: radiative driving and pulsa- 
tion driving. The arguments for these two possibilities 
are given in Neilson fc Lester (2008a), who also derive a 
model for pulsation-driving in Cepheids. It is not fea- 
sible to apply the pulsation-driving model to this set of 
data as we do not have knowledge of the pulsation ampli- 
tudes or masses to which the model is sensitive. However 
we can test whether a radiative-diiveTi stellar wind can 
match the predicted mass-loss r ates using the method of 
iCastor. Abbott, fc KleinI (fl975h . 

The calculation of the mass-loss rate for a radiative- 
driven stell ar wind is reviewed in iLamers fc Cassinellil 
dHH), and lNeilson fc Lested ^OOM) and will not be re- 
peated here. To conduct the calculation the mass, lumi- 
nosity, radius and effective temperature are needed; the 
radii are determined by the fitting, the effective tem- 
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Sieratures are given by the relation from iBeaulieu et al.l 
20011 ) and the luminosity is found from the radius 
and effective temperature. The mass is unknown so 
the radiative-driven mass-loss rates are found using a 
number of masse s via mass-luminosity relations from 
iBono et al.l (|2000f) . where L = M" in solar units. The 
mass-loss rates found from the observations are shown 
in the left panel of Figure [TOl plotted with the best- 
fit relations for radiative-driven mass-loss rates found 
with the following mass-luminosity relations, where n = 
4.4,4.7,5.0, and 6.0. The value of n = 4.4 represents 
the mass-luminosity relation from stellar evolution cal- 
culations, while n — 4.7 and 5.0 represent the mass- 
luminosity relations relating to mass found using pulsa- 
tion calculations, and n = 6.0 is used as an extreme case. 
An example of the values of the radiative-driven mass- 
loss rates is shown in the right panel of Figure (TU] for 
the case of n = 4.7. Radiative-driven mass-loss rates 
for other values of n > 4.7 will increase the rates and 
for n < 4.7 will decrease the rates by a roughly constant 
amount for each Cepheid. 

The radiative-driven mass-loss rates are significantly 
smaller than the mass-loss rates determined from the 
observations. At short periods of approximately 5 days, 
the radiative mass-loss rates are about 10"^ to 10^ times 
lower. However, at periods greater than 30 days the 
radiative-driven mass-loss rates are of similar order as 
the calculations. This implies that the mass-loss can- 
not be driven by radiative lines alone at short period; 
there must be another driving mechanism. This differs 
at longer period, but it should be noted that the mass- 
loss rates found from the observations are the minimum 
value based on the dust-to-gas ratio. This means that 
the predicted gas mass-loss rates from infrared observa- 
tions may be larger than the radiative-driYen mass-loss 
rates. 

As a further test of whether the mass loss is consistent 
with radiative driving we compute the circumstellar flux 
from dust created in a radiative-driven wind and added 
that to the predicted blackbody fluxes to compute in- 
frared PL relations. These relations are predictions of 
what would be observed if the mass loss is consistent 



with radiative driving. The fitted relations at the four 
wavelengths 3.6, 4.5, 5.8, and 8.0 fj.m have slopes and 
y-intercepts that are equivalent to the IR PL relations 
derived from the predicted stellar fluxes alone within the 
error of the fits. For instance the slope and y-interccpt of 
the 8.0 fim relation is -3.139 ±0.024 and 15.988 ±0.017, 
differing by only a few thousandths from the 8.0 nm 
PL relation determined from the predicted stellar fluxes 
alone. Radiative driving does not explain the significant 
infrared excess of 44 Cepheids that are explained by mass 
loss. 

The amount of mass loss from the LMC Cepheids do 
not agree with radiative-driving calculations, and the 
mass los s is more consistent with the pulsation-driven 
model of iNeilson fc Lesteil (|2008al ibl) if one considers the 
magnitude of the mass-loss rates and the amount of scat- 
ter, especially for the 44 Cepheids where the mass-loss 
model is unique. The remainder of the sample cannot 
be distinguished statistically from pulsation driving, ra- 
diative driving or no mass loss. However, if one invokes 
no mass loss or radiative driving then it is more difhcult 
to explain the behavior of the 44 Cepheids. Pulsation- 
driven mass loss predicts that mass loss is driven by 
shocks generated in a pulsating atmosphere and the 
shocks tend to be more efficient at hotter effective tem- 
peratures. This would explain the large range of mass- 
loss rates at similar pulsation periods. The pulsation- 
driven mass-loss rates tend to be orders of magnitude 
larger than radiative-driven mass-loss rates as shown for 
observations of Galactic Cepheids and theoretical models 
of Galactic, LMC and SMC Cepheids at shorter periods. 
We conclude that the mass-loss rates found in this work 
provide evidence for the model of pulsation-driven mass 
loss. 

6. DISCUSSION AND CONCLUSIONS 

In this article, we hypothesized that LMC Cepheids 
have significant infrared excess based on infrared ob- 
servations of the LMC from the SAGE survey, and 
that the infrared excess is caused by dust forming in 
a Cepheid wind at a large distance from the surface 
of the star. The idea was tested using OGLE II 
BVI observations along with IRAC observations of the 
Cepheids to best-fit radii and dust mass-loss rates of 
the LMC Cepheids. The effective temperatures are de- 
termined using a temp erature-color-period relation from 
IBeaulieu etldl (|2001l) . 

The predicted gas mass-loss rates are significant with 
an average about 10^^" to lO^^Mo/yr and may possibly 
be as high as lO'^A/o/yr, depending on the value of the 
dust-to-gas ratio. These mass-loss rates are not con- 
sistent with have radiative line-driving as the primary 
driving mechanism for the LMC Cepheids. The rates, 
instead, provide evidence for shocks and pulsation driv- 
ing the m ass loss when compared to the analytic model 
from Neils on fc Lesteil (2008a^b). 

The mass-loss model is compared to the fit of the ob- 
servations with radii alone, and it is shown that a model 
fitting only the radius may be rejected with 95% con- 
fidence relative to the mass-loss model for 44 Cepheids 
but this means that the remainder of the sample is con- 
sistent with no mass loss. Therefore the results and argu- 
ments in the work based on the mass-loss model should 
be regarded with caution. For almost every case the pre- 
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Fig. 10. — (Left) The comparison of the mass-loss rates found from the infrared observations with best— fit linear relations representing 
mass— loss rates found from radiative-driving calculations, and (Right) the radiative-driven mass— loss rates for the mass-luminosity relation 
where n = 4.7. 



dictions are limited by the uncertainties of the distance 
modulus and the IR pulsation amplitude. This compar- 
ison of the two models using the F-test is dependent 
on the mass-loss rate, and this shows we detect reliable 
infrared excess if the mass-loss rate is > 10~^°MQ/yr. 
Because of the dependence of the F-test on the mass-loss 
rate we explored what effect mass loss would have on the 
infrared PL relation. The mass-loss model would benefit 
from infrared observations over the period of pulsation 
to determine the mean brightness, which would decrease 
the uncertainty of the infrared brightnesses. 

The large mass-loss rates of LMC Cepheids may ex- 
plain the Cepheid mass discrepancy in the LMC. The 
pulsation masses tend to b e about 20% smaller than evo- 
lution masses in the LMC ()Kelleill2008l ). which translates 
to a difference of about IMq for lowest mass Cepheids 
up to a few solar masses for the most massive Cepheids. 
The mass-loss rates found in this work agree with the 
discrepancy for low mass Cepheids to an order of magni- 
tude because the evolutionary timescale for a Cepheid on 
its second crossing is of order ten million years. However 
the mass-loss rates are too small to be consistent with 
a 20% mass discrepancy for the more massive Cepheids. 
It should be noted that the mass discrepancy in LMC 
Cepheids is measured from Cepheids with periods less 
than 20d, which have evolu tionary masses from about 4 
to 7Mq (|Bono et al.ll2000[) . The mass discrepancy has 
not been measured for more massive LMC Cepheids. 

It has also been found that mass loss affects the in- 
frared PL relations. Using the predicted stellar luminosi- 
ties we constructed new infrared PL relations that do not 
have infrared excess. These relations differ from those 
determined bv lNgeow fc Kanburl (|2008( ) with differences 
in th e zero point and t he slo pe. The IR PL relations 
from iNgeow fc Kanbuil (2008) have slopes that become 
smaller at longer wavelength inconsistent with the argu- 
ment that the slope of the PL relation should approach 
a constant maximum value at longer wavelength bas e 
on the Period-Radius relation (jFreedman et al.l [2008D . 



The slopes in this work are all about —3.15 with a small 
amount of steepening at longer wavelength. This would 
imply a constant slope near that value which is also in- 
consistent with the slope derived from the PR relation. 

Using the F-test, there is evidence for non-linearity 
in the relations similar to the non-linear structure found 
in optical PL relations. Mass loss acts to linearize the 
PL relation at 3.6,4.5, and 5.8 fim, while at 8.0 fim the 
PL relation is non-linear with a slope that is shallower 
at P < lOd than for P > lOd which implies the infrared 
excess is becoming more important at longer wavelength. 
Mass loss may also explain why the K-band PL relation 
IS marg inally hnear lNgeow fc Kanburl ()2006[ ). 

The resulting effect that mass loss has on infrared ob- 
servations of Cepheids implies serious consequences for 
infrared Period-Luminosity relations if they are to be 
used for high precision astrophysics. One of the rea- 
sons for using infrared PL relations is that they are 
less sensitive to metallicity than optical relations and 
hence do not need to be corrected for each galaxy 
(jSasselov et al.l I1997D . The metallicity correction is a 
significant so urce of uncertainty in studies of the Hub- 
ble Constant (jFreedman et al.ll2001[ ) and an infrared PL 
relation that avoids this uncertainty would be a power- 
ful tool. However, we have shown that mass loss affects 
the scatter and the structure of the PL relation. The 
scatter increases the uncertainty of any distance deter- 
mination, but more importantly the fractional amount 
of dust generated in a wind depends on metallicity. The 
amount of mass los s may also depend on metallicity, as 
suggested by Neil son fc Lesteil ()2008bl ). These two is- 
sues imply the Period-Luminosity relation depends on 
metallicity at infrared wavelengths as well as at optical 
wavelengths though to what extent is currently unknown. 
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